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Abstract 

We examine whether an appreciable value of the lepton flavor mixing angle 
613 at the electroweak scale Aew can be radiatively generated from 613 = 0° 
at the GUT scale Aqut- It is found that the renormalization-group running 
and seesaw threshold effects may lead to #13 ~ 5° at low energies for two 
simple large-maximal mixing patterns of the MNS matrix in the minimal 
supersymmetric standard model. If #12 is sufficiently large at Agut> it will 
be possible to radiatively produce #13 ~ 5° at Aew both in the standard 
model and in its supersymmetric extensions. The mass spectrum of three 
heavy right-handed Majorana neutrinos and the cosmological baryon number 
asymmetry via leptogenesis are also calculated. 
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I. INTRODUCTION 



The recent solar [1], atmospheric [2], reactor (KamLAND [3] and CHOOZ [4]) and accel- 
erator (K2K [5]) neutrino oscillation experiments have provided us with very robust evidence 
that neutrinos are massive particles and their mixing involves two large angles (0i 2 ~ 33° 
and #23 ~ 45°) and one small angle (0i3 < 13°). How small 0i3 is remains an open question, 
but a global analysis of the presently available neutrino oscillation data [6] indicates that #13 
is most likely to lie in the range 4° < #13 < 6°. One important target of the future neutrino 
experiments is just to measure #13 [7]. 

The smallness of 6> 13 requires a good theoretical reason, which might simultaneously 
account for the largeness of #12 and #23- If #13 = 0° held, there should exist a kind of 
new flavor symmetry which forbids flavor mixing between the first and third lepton families. 
While such a new symmetry is unlikely to exist at or below the electroweak scale (Aew ~ 10 2 
GeV), it might show up at a superhigh scale - e.g., the scale of grand unified theories 
(Aqut ~ 10 16 GeV). Then a natural way to break this flavor symmetry and obtain 13 7^ 0° 
in a specific model is to run relevant parameters of the model from Aqut to Aew by making 
use of the renormalization group equations (RGEs) [8] and taking account of the seesaw 
threshold effects [9], either in the standard model (SM) or in the minimal supersymmetric 
standard model (MSSM). Antusch et al have recently presented two simple examples (one 
with 6»i2 = #23 = 45° and 13 = 0° at A GUT [10], and the other with 12 = 13 = 0° and 
023 = 45° at Aqut [11]) to radiatively generate 9 12 ~ 33° and #13 7^ 0°, but their primary 
interest is in 9 12 and their results for 0i3 are far below the best-fit values of 0i3 obtained 
from the global analysis [6]. Although the low-scale output of 0i3 is somehow adjustable by 
scanning the parameter space of a given model at the GUT scale, we find that it is highly 
nontrivial to obtain 0i3(A EW ) ~ 5° from 0i3(A GUT ) = 0° and fit all experimental data of 
neutrino oscillations in the meantime. 

The main purpose of this paper is to examine whether an appreciable magnitude of 
013 can be radiatively generated, from Aqut to Aew? through the seesaw thresholds. We 
shall consider four instructive patterns of the Maki-Nakagawa-Sakata (MNS) lepton mixing 
matrix Vmns at A G ut as typical examples: 
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where P$ = Diag{l, 1, e lS } with 5 being a CP-violating phase in the standard parametrization 
of Vmns [12] 1 . While patterns (A) and (B) with CP conservation (i.e., 5 = 0°) have been 
discussed in Refs. [10] and [11], we shall show that S can actually play an important role in 
the radiative generation of 6> 13 . Patterns (C) [13] and (D) [14] are phenomenologically favored 
to account for current experimental data of solar and atmospheric neutrino oscillations. We 
find that the RGE running and seesaw threshold effects may allow us to obtain #13 ~ 5° at 
low energies from both patterns (C) and (D) in the MSSM. If 9 12 ~ 60° is taken at Aqut, it 
will be possible to radiatively produce 9± 3 ~ 5° at A EW both in the SM and in the MSSM. As 
a by-product, the mass spectrum of three heavy Majorana neutrinos and the cosmo logical 
baryon number asymmetry via leptogenesis are also calculated. 

II. RGE RUNNING AND THRESHOLD EFFECTS 

Let us make a simple modification of the SM by introducing three heavy right-handed 
neutrinos iVj (for i = 1,2,3) and keeping the Lagrangian of electroweak interactions invariant 
under SU(2) L x U(l)y gauge transformation. In this case, the Lagrangian relevant for lepton 
masses can be written as 

-A epton = TlXie K H + l L Y u U R H C + ^MrI/r + h.c. , (2) 

where /l denotes the left-handed lepton doublet; e R and v K stand respectively for the right- 
handed charged lepton and Majorana neutrino singlets; and H is the Higgs-boson weak 
isodoublet (with H c = ia 2 H*). If the MSSM is taken into account, one may similarly write 
out the Lagrangian relevant for lepton masses: 

-lepton = ^e R #! + l L Y u u R H 2 + ^M R v R + h.c. , (3) 

where Hi and H 2 (with hypercharges ±1/2) are the MSSM Higgs doublets. To obtain the 
effective (left-handed) neutrino mass matrix, a common approach is to integrate M R out 
of the full theory. This corresponds to a replacement of the last two terms in £i cpt on or 
^Lpton by a dimension-5 operator, whose coupling matrix takes the well-known seesaw form 
k = —Y v M R l Yj [15]. However, the threshold effects in such a naive treatment have to be 
taken into account, because the mass eigenvalues of Ni may have a strong hierarchy (for 
example, M 3 > M 2 > M x ). 

We take Aepton or C[ cpton for granted at the GUT scale, where the Yukawa interactions 
of quarks and Higgs bosons with the coupling matrices Y u (up) and Y d (down) can similarly 
be written out. To evolve the lepton mixing parameters from Aqut to Aew m a generic 
seesaw model 2 , one has to make use of a series of effective theories which are obtained by 



1 For simplicity and illustration, we only take account of a single CP-violating phase in this work. 

2 We assume the supersymmetry breaking scale Asusy to be close to the electroweak scale Aew> 
just for the sake of simplicity. Even if Asusy/Aew ~ 10 holds, the relevant RGE running effects 
between these two scales are negligibly small for the physics under consideration. 
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integrating out the heavy right-handed singlets JVj step by step at their mass thresholds. 
The derivation of the one-loop RGEs and the method for dealing with the effective theories 
have been presented in Ref. [9] in an elegant way. Here we summarize a few essential steps 
to be taken in treating the seesaw threshold effects. 

(a) We use the one- loop RGEs to run Y v and M R from A G ut to the heaviest right-handed 
neutrino mass scale M 3 . A proper unitary transformation of the right-handed neutrino fields 
allows us to diagonalize M R at M 3 - namely, U^M R U^ = Diag{Mi, M 2 , M 3 }. Then Y v is 
transformed into Y U U^. The effective neutrino coupling matrix can be obtained by 
integrating out M 3 . In this case, we denote 3 



Y u(3) = Y U U, 
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Y m = Y„Ur I I , M R(3) = ( ^ ^ ) (4) 



and get the tree-level matching relation k, 3 n = — Y U ^M% 1 Yj^, where all variables have been 
set to the scale \x = M 3 . 

(b) We further run Y^ 3 ) , M R ( 3 ) and from M 3 down to the intermediate right-handed 
neutrino mass scale M 2 . Because the RGE running effect may spoil the diagonal feature 
of M R ( 3 ), a re-diagonalization of M R ( 3 ) at M 2 should be done by means of a 2 x 2 unitary 
transformation matrix C/ R . Integrating out M 2 , we arrive at 

Yu(2) = n ( s)tfc(J) , t( 2) = n ( s)t%(J) , M R(2) = (Mi) (5) 

and the tree- level matching condition = — Y^M^Yj^, where all variables have 
been set to the scale \i = M 2 . 

(c) We follow a similar way to run Y u ( 2 ), ^r(2) and from M 2 down to the lightest 
right-handed neutrino mass scale M x . As M R ( 2 ) is actually a 1 x 1 mass matrix, it does not 
need to be re-diagonalized at M\. Integrating out Mi, we obtain 

k = = « (2) - Y v{1) M^Yj {1) , (6) 

where = ^(2) holds, and all variables have been set to the scale /x = Mi. 

(d) Finally, we run k from M x down to the electroweak scale A EW . The one- loop RGE 
governing the evolution of k is given by [8] 



167T 2 — — = an + C 
at 



(y 1 Y^)k + k(y 1 Y^ T ] , (7) 



where t = ln(/x/Mi) with \i being the renormalization scale. We have C = —1.5, a ~ 
-3#f + 6/ t 2 + A in the SM and C = 1, a « -1.2#? - 6#f + 6/ 2 in the MSSM [16], where 
g 12 denote the gauge couplings, f t denotes the top-quark Yukawa coupling, and A denotes 



3 It should be noted that our notations (in particular, k = —Y^M^Yj) are somehow different 
from those of Ref. [9], where the seesaw formula k = 2Yj M^ l Y u has been adopted. 
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the Higgs self-coupling in the SM. After spontaneous gauge symmetry breaking, we arrive 
at the fermion mass matrices M v = v 2 k, Mi = vYi, M u = vY n and M<j = vY d in the SM; and 
M v = v 2 k sin 2 0, Mi = vYi cos (3, M u = vY u sin (3 and M d = vY d cos (3 in the MSSM, where 
v ks 174 GeV stands for the vacuum expectation value of the neutral Higgs field in the SM, 
and tan (3 represents the ratio of two vacuum expectation values in the MSSM. 



III. INITIAL CONDITIONS AND ASSUMPTIONS 



The lepton (or quark) flavor mixing matrix Vmns ( or Vckm) arises from the mismatch 
between the diagonalizations of Y\ (or Y u ) and k (or Yd). Without loss of generality, we 
arrange Y% and Y u to be diagonal, real and positive at A GU t; i-e., 
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where Qi = ^2 = v in the SM, and Q\ = v sin (3 and ^2 
flavor basis, Y& and Y v can generally be expressed as 



1 



■v ( 



CKM 



m d 






















m 



U d , Y v 



v cos (3 in the MSSM. In this 



(9) 




where n, r 2 and y u are three real and positive dimensionless parameters characterizing the 
eigenvalues of Y v \ and V u , U v and are three 3x3 unitary matrices. It is obvious that 
the complex phases of V v and U v can be re-arranged, such that the former contains a single 
irremovable CP-violating phase as Vqkm does. Furthermore, one may redefine the relevant 
right-handed fields of quarks and neutrinos to rotate away from Y d and U v from Y v . Such 
a transformation of Y u is equivalent to absorbing U v into the Majorana mass matrix M R , 
which is not required to be diagonal at A GU t- Once and U v are rejected, we are only 
left with seven unknown parameters in the above Yukawa coupling matrices (namely, three 
eigenvalues of Y v and four mixing parameters of V v ). 

To fix the pattern of Mr at Aqut, we extrapolate the effective neutrino coupling matrix 
k and the lepton flavor mixing matrix Vmns up to the GUT scale: 
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where mj (for i = 1, 2, 3) denote the physical masses of three light neutrinos. Then M R can 
be determined from the inverted seesaw formula M R = -Yj k~ x Y v [19]. Note that Vmns 
consists of three mixing angles and three CP-violating phases. For the sake of simplicity, 
here we only take account of the "Dirac-like" phase of Vmns- Then the unitary matrices 
Vckm, Vmns and V v may universally be parametrized as 
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where Cij = cos% and Sij = 9ij. In the limit of #13 = 0°, the complex phase of Vmns actually 
serves as a "Majorana-like" phase and may significantly affect the RGE running behaviors 
of neutrino masses and lepton flavor mixing angles [17]. Four typical patterns of Vmns 5 as 
already listed in Eq. (1), will be taken in our following investigation. 

To be specific, we only pay attention to the normal mass hierarchy of three light neutrinos 
(i.e., m 3 > m 2 > mi). Then we arrive at m 2 = \Jm\ + Am^ m and m 3 = \jm\ + Aml tm , 
where Am^ and Am^ tm denote the mass-squared differences of solar and atmospheric 
neutrino oscillations. Since the low-scale values of (m u ,m c ,m t ), (m<i,m s ,mb), (m^m^r/v) 
and (Am^jAm^) are all known, we just have a single unknown mass parameter (mi). 
On the other hand, four parameters of Vckm are also known at low energies [18]. Given 
a special pattern of Vmns m Eq. (1), only its CP-violating phase is not fixed. In short, 
we are totally left with nine free parameters at A GUT : the lightest neutrino mass mi, three 
eigenvalues of Y v , four mixing parameters of V v and the CP-violating phase of Vmns- We 
should also specify the value of the Higgs mass m H in the SM and that of tan (3 in the MSSM 
when numerically solving the relevant RGEs. 

Although y v , r\ and r 2 are arbitrary parameters, they are expected to be of or below 
Oil). The unknown rotation and phase angles of V u and Vmns are allowed to take possible 
values in the range between and 2n. As for those parameters whose sizes are known at 
low energies, one may properly adjust their initial values at A GUT and run the RGEs to 
reproduce their low-scale values within reasonable error bars. We are then able to fix the 
parameter space by fitting all relevant low-scale data. Such a phenomenological approach 
will allow us to examine whether 6*13 (Aew) ~ 5° can be generated from #13 (Agut) = 0° for 
Vmns- Naively, we speculate that an appreciable RGE enhancement of 13 may take place 
either between the scales Agut and M 1 or between the scales Mi and Aew- m either case, 
the masses of three light neutrinos are required to be nearly degenerate. As the seesaw 
threshold effect can significantly affect the RGE running behaviors in most cases [9-11,20], 
it is more likely that the dominant RGE enhancement of #13 occurs between Agut and Mi. 
We shall demonstrate this observation in our subsequent numerical calculations. 

IV. NUMERICAL EXAMPLES AND DISCUSSIONS 

We carry out a numerical analysis of the RGE running and seesaw threshold effects by 
following the strategies outlined above and taking four typical patterns of Vmns at the GUT 
scale. Our results are summarized in Table 1 and Figs. 1 and 2. Some comments and 
discussions are in order. 

(1) We have carefully examined the sensitivity of Vmns (Aew) to the value of every free 
parameter at Agut- We find that the inputs of mi, y v , #^> NS , 5mnsj #13 °f V v and Am^ 
appear to be very important in adjusting the output of #^ NS at A E Wj while the other 
parameters mainly play a role in fine-tuning the results. In particular, the effects of r 1 , r 2 , 
9 12 of V v , 923 °f Vv and 5 of V v are insignificant, because the dominant RGE enhancement of 
#^| NS takes place above the heaviest right-handed neutrino mass scale M3. Hence we have 
simply fixed 6*12 = #23 = S = 0° for V u at Agut- Note that the initial values of <5mns and 
#13 of V u are important in specifying the running tendency of two neutrino mass-squared 
differences and three mixing angles of Vmns 5 and they have to be sufficiently large (as shown 
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in Table 1) in order to generate an appreciable magnitude of 9\™ s at low energies. To 
be more explicit, the evolution of 6^ NS strongly relies on <5mns m the fitting, while the 
enhancement of #^ NS (from zero at A GU t to a few degrees at A EW ) requires a big input 
value for 9 13 of V u (~ 45°, for example). It is also worth remarking that the parameter space 
obtained here should by no means be unique; or rather, it mainly serves for illustration. To 
exhaustively explore the allowed ranges of all relevant parameters is a quite lengthy work 
and will be presented elsewhere [21]. 

(2) Restricting ourselves to the typical parameter space illustrated in Table 1, we find 
that it is difficult to produce 9^ NS ~ 5° at A EW in the SM. Our results yield 9^ NS ~ 3° 
for pattern (A) and 6^ NS ~ 1.5° for patterns (B), (C) and (D) at low energies. In the 
MSSM, however, a more strong RGE enhancement of 0^ NS becomes possible. For example, 
we obtain 6>^ NS ~ 3° for pattern (B) with tan/3 ~ 10 and 6»^ NS ~ 5° for patterns (C) 
and (D) with tan/3 ~ 20 at the electroweak scale. We think that the numerical results in 
the supersymmetric case are encouraging for model building, because a new kind of flavor 
symmetry at A GUT might naturally lead to the bi-maximal mixing pattern (B) or the large- 
maximal mixing patterns (C) and (D). 

(3) The running behavior of #^ NS deserves some more remarks. As already noticed in 
Refs. [8-1 1] , the evolution of #^, NS from A EW to A GU t ( or vice versa) depends sensitively on 
how big its initial value is and whether three light neutrinos have a near mass degeneracy. 
Figs. 1 and 2 illustrate that #^ NS (A EW ) ~ 33° can radiatively be obtained, either from 
#^ ns (Agut) = 0° in pattern (A) or from ^ NS (A GUT ) ~ (30° -45°) in patterns (B), (C) and 
(D). For pattern (A) or (B), the RGE running of 6»^ NS in the MSSM is quite similar to that 
in the SM. The reason for this similarity is simply that the magnitudes of mi (SM) and mi 
(MSSM) at A EW are comparable (~ 0.05 eV) and the value of tan/3 in the MSSM case is 
mild (~ 10). When patterns (C) and (D) are concerned, however, the running behavior of 
0^ NS in the MSSM is more violent than that in the SM. The reason for such a remarkable 
difference is two-fold: first, mi(A E w) ~ 0.14 eV and mi(A G ijr) ~ 0.20 eV in the MSSM case 
imply that the masses of three light neutrinos are nearly degenerate in the entire running 
course - this near mass degeneracy can significantly affect the behavior of 9f™ s [8]j second, 
the large value of tan (3 (~ 20) plays an important role in enhancing the RGE evolution of 
light neutrino masses and flavor mixing angles (e.g., 0^ NS oc (1 + tan 2 /?) [8]). 

(4) We stress that 9^ NS ~ 5° at A E w can be radiatively generated from 6>J| NS = 0° at 
A G ut even in the SM, only if we go beyond the four simple patterns considered above. A key 
point is to assume 6^ NS to be large enough at the GUT scale. Such an example is presented 
in Fig. 3, where 9™ NS = 67° (SM) versus 9^ NS = 63° (MSSM) has been taken. One may see 
that the dominant RGE suppression of #f 2 NS occurs from A GU t to M 3 , and the dominant 
RGE enhancement of #J| NS takes place in the same region. Why #J^ NS > 45° holds at A G ut 
is certainly a big question. Putting aside this question, we remark that #^ NS ~ 33° at low 
energies can in principle be produced from an arbitrary value of #^ NS at the GUT scale 
via the seesaw threshold effects. Furthermore, the radiative generation of 9^ NS is highly 
sensitive to the initial condition of #j^ NS . We observe that the running of #^ NS is rather 
stable, unlike 6»^ NS and 6»^ NS . Note that the CP-violating phase <5mns is also stable against 
radiative corrections from A G ut to A E w- This conclusion is true for both the example under 
consideration and the four patterns discussed above. 

(5) A by-product of our analysis is the determination of three heavy right-handed neu- 



7 



trino masses, as shown in Table 1. We see that they have a clear normal hierarchy. It is 
then possible to calculate the cosmological baryon number asymmetry via leptogenesis [22]. 
Instead of describing the technical details of leptogenesis [23], we only give a brief summary 
of its essential points in the following. Lepton number violation induced by the third term 
of Aepton in Eq. (1) or of C[ epton in Eq. (2) allows decays of three heavy Majorana neutrinos 
Ni to happen: Ni — > I + h and Ni — > I + h c , where h = H in the SM or h = H 2 in the 
MSSM. Because each decay mode occurs at both tree and one-loop levels (via the self-energy 
and vertex corrections), the interference between these two decay amplitudes may result in 
a CP- violating asymmetry between —> I + h and its CP-conjugated process. If the 
interactions of Ni are in thermal equilibrium when A3 and N 2 decay, the asymmetries e 3 
and e 2 can be erased before N\ decays. Then only the asymmetry E\ produced by the out- 
of-equilibrium decay of N\ survives. This CP- violating asymmetry may lead to a net lepton 
number asymmetry Yl oc e±, and the latter is eventually converted into a net baryon number 
asymmetry Y B via the nonperturbative sphaleron processes [24]: Y B ps — 0.55Yl in the SM 
or Y B — 0.53F L in the MSSM. We follow these steps to evaluate Y B for four patterns listed 
in Table 1. It turns out that only pattern (B) can yield an appreciable result, which lies in 
the generous range 0.7 x 1(T 10 < Y B < 1.0 x 10~ 10 drawn from the recent WMAP obser- 
vational data [25]. Compared with pattern (B), patterns (A), (C) and (D) have relatively 
small values of M 1 (~ 10 9 GeV) and M 1 /M 2 (~ 3 x 10~ 2 ). The outputs of Y B in these three 
patterns are therefore suppressed and below the observational result. If the example shown 
in Fig. 3 is taken into account, we may obtain Y B ps 7.7 x 10~ n (SM) or Y B ps 8.0 x 10 -11 
(MSSM), which is compatible with the WMAP data. 

Finally, it is worthwhile to point out that we have also taken a look at the "democratic" 
neutrino mixing pattern [26] 
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at the GUT scale and examined the possibility to radiatively produce # 13 ~ 5° at the 
electroweak scale and to simultaneously fit all experimental data of neutrino oscillations. 
We find that it is extremely difficult, if not impossible, to obtain #23 ~ 45° at low energies. 
The reason is simply that the initial value of 23 (ps 54.7°) is not close to 45° and the RGE 
evolution cannot significantly change the magnitude of this mixing angle. Therefore, we 
argue that radiative corrections might not be a very natural way to break the lepton flavor 
democracy. The latter could be explicitly broken in some other possible ways at or below 
the GUT scale. 



V. SUMMARY 

We have conjectured that the lepton flavor mixing angle 9 13 might be vanishing at the 
GUT scale due to the existence of a kind of new flavor symmetry. Starting from this point 
of view, we have examined whether an appreciable value of #13 at low energies can be radia- 
tively generated via the RGE running and seesaw threshold effects. Four simple but typical 
patterns of the lepton flavor mixing matrices have been considered as the initial conditions 
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at Agut- It is found that the dominant RGE enhancement of #13 takes place from Aqut 
to the heaviest right-handed neutrino mass scale. For two large-maximal mixing patterns, 
#i3(Aew) ~ 5° can be radiatively produced in the MSSM. We have also demonstrated that 
it is possible to obtain 6> 13 ~ 5° at A EW from 9 13 = 0° at A GUT in the SM, if the initial value 
of #12 is large enough. As a useful by-product, the mass spectrum of three heavy Majorana 
neutrinos is determined and the cosmological baryon number asymmetry via leptogenesis is 
calculated. 

Although the numerical examples presented in this paper are mainly for the purpose of 
illustration, they are quite suggestive for model building. Of course, only the future neutrino 
experiments can tell us how small 6> 13 is. But we believe that the radiative generation of 
#13 from high to low energies is an interesting theoretical approach towards understanding 
the smallness of #13, and it might indicate some useful hints about the underlying flavor 
symmetry which is associated with the dynamics of lepton flavor mixing and CP violation. 

This work was supported in part by the National Nature Science Foundation of China. 
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TABLES 



TABLE I. Numerical examples for radiative generation of #13 via the RGE evolution and seesaw 
threshold effects. The cosmological baryon number asymmetry Yq is also computed. Note that we 
have taken 9 12 = #23 = 8 = 0° for V v at Aqut- 
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FIG. 1. The RGE evolution of three lepton mixing angles between Aew an d Agtjt m the SM, 
where the initial values of relevant parameters are listed in Table 1. 
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FIG. 2. The RGE evolution of three lepton mixing angles between Aew and Aqut in the 
MSSM, where the initial values of relevant parameters are listed in Table 1. 
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FIG. 3. Running of three lepton mixing angles, where mi = 0.15 eV, Am s 2 
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7.6 xlO- 3 eV 2 , y„ = 0.9, n = 1/41, r 2 = 1/17, {0 12 , 23 ,0i 3 , <5} yMNS = {67°, 45°, 0°, 94°}, 
{#12, 6*23, 6»i 3 , 5}^ = {0°, 0°, 45°, 0°} and m H = 120 GeV have typically been input at Aqut in the 
SM; and m-y = 0.12 eV, Am 2 un = 3.9 x 10" 4 eV 2 , Am 2 tm = 5.4 x 10~ 3 eV 2 , y v = 0.8, n = 1/41, 
r 2 = 1/15, {0 12 , 23 ,0i 3 , 5} Vmns = {63°, 45°, 0°, 100°}, {6 12 ,9 23 ,9 13 ,5} v „ = {0°, 0°, 45°, 0°} and 
tan (3 = 10 have typically been input at Agut in the MSSM. 
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